On the Rational Terms of the one-loop 

amplitudes 



Giovanni OssoW 1 ) , Costas G. Papadopoulos^ , Roberto Pittau^ 2,3 ^ 

(^Institute of Nuclear Physics, NCSR Demokritos, 15310 Athens, Greece 

^Departamento de Fisica Teorica y del Cosmos 
Centro Andaluz de Fisica de Particulas Elementales (CAFPE) 
Universidad de Granada, E-18071 Granada, Spain 

^Dipartimento di Fisica Teorica 
Univ. di Torino and INFN sez. di Torino 
V. P. Giuria 1, 1-10125 Torino, Italy 



Abstract 

The various sources of Rational Terms contributing to the one-loop amplitudes 
are critically discussed. We show that the terms originating from the generic (n — 
4)— dimensional structure of the numerator of the one-loop amplitude can be derived 
by using appropriate Feynman rules within a tree-like computation. For the terms 
that originate from the reduction of the 4— dimensional part of the numerator, we 
present two different strategies and explicit algorithms to compute them. 



1 



Contents 

1 Introduction 

2 The origin of R 2 H 

3 The OPP method and the origin of Ri Q 

4 The n-dimensional cuttings to compute Ri 

5 Conclusions 
A The g 2 dependence of the OPP coefficients 



2 



1 Introduction 



In the last few years, a big effort has been devoted by several authors to the problem of 
computing one-loop amplitudes efficiently [1]. Besides Standard techniques, where ten- 
sor reduction/computation is performed, numerically or analytically, new developments 
emerged, originally inspired by unitarity arguments (the so called unitary and generalized 
unitarity methods) [2] - [5], in which the tensor reduction/computation is substituted by 
the problem of determining the coefficients of the contributing scalar one-loop functions. 
This possibility relies on the fact that the basis of one-loop integrals is known in terms 
of Boxes, Triangles, Bubbles and (in massive theories) Tadpoles, in such a way that, 
schematically, one can write a Master Equation for any one-loop amplitude M. such as: 

M. = dt Boxj + c i Triangle j + 6j Bubble j + a.; Tadpole^ + R , (1) 

i i i i 

where di, Ci, bi and a. t are the coefficients to be determined. 

Nevertheless, in practice, only the part of the amplitude proportional to the one- 
loop scalar functions can be obtained straightforwardly in the unitary cut method. The 
remaining Rational Terms, RTs, (R in the above equation) should be reconstructed by 
other means, either by direct computation [6,7] or by boostrapping methods [8]. In [9] the 
RTs are obtained by explicitly computing the amplitude at different integer value of the 
space-time dimensions. On the other hand, in another recently proposed method, QPP [10], 
a class of terms contributing to R can be naturally derived in the same framework used 
to determine all the other coefficients. 

In this paper, we critically analyze the various sources of RTs appearing in one-loop 
amplitudes, by classifying them in two categories: R = Ri + R 2 , also presenting a few 
computational methods. In the next section, we investigate the origin of R 2 and develop 
a practical computational strategy. In Section [3[ after a brief recall of the OPP method, 
we present a way to compute Ri, strictly connected to the OPP framework. In Section BJ 
we describe yet another method that relies on cuts in n-dimensions. This second method 
is proven more suitable for numerical applications within the OPP algorithm and may 
also be applied in a more general framework. Finally, in the last section, we outline our 
conclusions. 



2 The origin of R2 

Our starting point is the general expression for the integrand of a generic m-point one-loop 
(sub-) amplitude [10] 

A ® = n n Ni<l) n ' A = (q + Pif -m 2 , Po ? . (2) 

In the previous equation, dimensional regularization is assumed, so that we use a bar 
to denote objects living in n = 4 + e dimensions. Furthermore q 2 = q 2 + g 2 , where q 2 
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is e-dimensional and (q • q) = 0. The numerator function N(q) can be also split into a 
4-dimensional plus a e-dimensional part 



N(q) = N(q)+N(q\q,e). (3) 

N(q) is 4-dimensional (and will be discussed in the next section) while N(q 2 ,q,e) gives 
rise to the RTs of kind R 2 , defined as 

To investigate the explicit form of N(q 2 ,q,e) it is important to understand better the 
separation in Eq. ((3]). From a given integrand A(q) this is obtained by splitting, in the 
numerator function, the n-dimensional integration momentum q, the n-dimensional 7 
matrices 7^ and the n-dimensional metric tensor into a 4-dimensional component plus 
remaining pieces: 

q = q + q, 

gP> = g>" + gi*. (5) 

Notice that, when a n-dimensional index is contracted with a 4-dimensional (observable) 
vector v^, the 4-dimensional part is automatically selected. For example 

q ■ v = q ■ v and f = f . (6) 

A practical way to compute R2 is determining, once for all, tree-level like Feynman 
Rules for the theory at hand by calculating, with the help Eq. ([5]) , the R2 part coming from 
one-particle irreducible amplitudes up to four external legs. The fact that four external 
legs are enough is guaranteed by the ultraviolet nature of the RTs, proven in [6]. 

As an illustrative example, we derive the complete set of the Feynman Rules needed 
in QED. Along such a line R 2 can be straightforwardly computed in any theory. We start 
from the one-loop ^e + e~ amplitude in Fig. CD The numerator can be written as follows 

N(q) = e 3 {75 (#1 + m e ) 7m (#2 + m e ) ^} 
= e 3 + m e )7 M ($ 2 + m e ) 7 /3 

- e ($1 - m e )7 M ($ 2 - m e ) + eq 2 7^ - q 2 7/37 M 7 /3 } , (7) 

where all e-dimensional 7-algebra has been explicitly worked out in order to get the 
desired splitting. The first term in the l.h.s of Eq. ([7|) is N(q), while the sum of the 
remaining three define N(q 2 , q, e) for the case at hand. By inserting N(q 2 , q, e) in Eq. (PJ 



1 e-dimensional 7 matrices freely anti-commute with four-dimensional ones: {7^,71/} = 0. 
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Qi = q + v\ = Qi + q 
Q2 = q + P2 = Q2 + q 



Do 

Di = {q + Pi) 2 
D 2 = {q + P2? 



Figure 1: QED •ye + e diagram in n dimensions, 
and using the fact that 



^5ffe - + (8) 



gives 



ie 3 



K 2 = -— 1 , + 0{e) 1 (9) 



57T 



that can be used to define the following effective vertex: 




Figure 2: QED / ye + e effective vertex contributing to R 2 . 

With analogous techniques, taking also into account the integrals given in Eq. ( TTTll . 
one gets all the remaining QED effective vertices given in Fig. [3l 

To summarize, the problem of computing R2 is reduced to a tree level calculation and 
we consider it fully solved. The Ri part is, instead, deeply connected to the structure of the 
one-loop amplitude, as we shall see in the next section. It is worthwhile to mention that 
only the full R = Ri + R2 constitutes a physical gauge-invariant quantity in dimensional 
regularization. 



3 The OPP method and the origin of Ri 

The OPP reduction algorithm provides a useful framework to understand the origin of the 
RTs of kind Ri. The starting point of the method is an expansion of N(q) in terms of 
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4-dimensional denominators A = (q + pi) 2 — m 2 

m—1 m—1 

N (<l) = [d^oH^is) + d{q; ioi^h)] Yl D ' 

m—1 m—1 

+ [ c {ivhh) + c(q; i iii 2 )} J[ A 

«0<«1<*2 i^J0i*li*2 
m—1 m—1 

+ [6(io*i) + 6(g;*o*i)] II A 

io<i\ i^io,ii 
m—1 m—1 

+ t a (^o) + a(q; i )) J[ A 

io j^i 
m—1 

+ II A • (10) 

i 

Inserted back in Eq. (j2j), this expression simply states the multi-pole nature of any m- 
point one-loop amplitude. The fact that only terms up to 4 poles appear is due to the 
fact that m-point scalar loop functions with m > 4 are always expressible in terms of 
boxes up to contributions 0(e). The last term with no poles, P(q), has been inserted 
for generality, but is zero in practical calculations where m-point amplitudes behave such 
as N(Xq) — > X m when A — » oo. The coefficients of the poles can be further split in two 
pieces. A piece that still depend on q (the terms d, c, b, a), that vanishes upon integration 
due to Lorentz invariance, and a piece that do not depend on q (the terms d, c, 6, a). Such 
a separation is always possible, as shown in [10], and, with this choice, the latter set of 
coefficients is therefore immediately interpretable as the ensemble of the coefficients of all 
possible 4, 3, 2, 1-point one-loop functions contributing to the amplitude. 

Once Eq. (fTOl) is established, the task of computing the one-loop amplitude is then 
reduced, in the OPP method, to the algebraical problem of fitting the coefficients d, c, b, a 
by evaluating the function N(q) a sufficient number of times, at different values of q, and 
then inverting the system. Notice that this can be performed at the amplitude level and 
that one does not need to repeat the work for all Feynman diagrams, provided their sum 



ie- 
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g^(2mi-p 2 /3) 



P 



a 



le 




16tt 2 



-p + 2m e 



Figure 3: QED 77, ee and 7777 effective vertices contributing to R2. 
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is known. 

The OPP expansion is written in terms of 4-dimensional denominators. On the other 
hand, n-dimensional denominators A appear in Eq. j2j), that differ by an amount q 2 from 
their 4-dimensional counterparts 

A = A + f (11) 

The result of this is a mismatch in the cancellation of the n-dimensional denominators 
of Eq. (J2j) with the 4-dimensional ones of Eq. ( fit)]) (the OPP expansion), that originates a 
Rational Part. In fact, by inserting Eq. (jTTT) into Eq. (fTOl) . one can rewrite it in terms of 
n-dimensional denominators (therefore restoring the exact cancellation), but at the price 
of adding an extra piece f(q 2 , q). The RTs of kind Ri are defined as 

^ikl^ w^h - (12) 

The explicit form of the function f(q 2 ,q) can be easily and explicitly obtained, in 
the framework of the OPP method, by rewriting any denominator appearing in Eq. (J2J) as 
follows 

k-w mth ^('-e)' (13) 



This results in 



A(q) = N{q) Z Z 1 ■ ■ ■ Z m _ x + K 2 , (14) 

UqU\ ■ ■ ■ L>m-l 



where TZ 2 is the integrand function introduced in Eq. (jlj). Then, by inserting Eq. (jTOI ) in 
Eq. (JHJ), one obtains 

t/_x ^ d^pi^h) + d{q;i iii 2 h) "j-r 1 = 

~ ^ D-D- D D H i 

^ c(w 2 ) + c(q; w' 2 ) ^ ^ 
+ 2- /; /; /; 11 Z ' 

«0<*1<*2 *° 11 * 2 17^01*11*2 

^ &(¥l)+%i¥l) Tr* * 

+ 5-5: 11 z i 

io<ii l o *i i^jo,ii 

+ ^ M*o) + a(g;*o) jjU 

m— 1 

+ P(g) \{2 l + n 2 . (15) 

Ri is then produced, after integrating over <i n g, by the q 2 dependence coming from the 
various 2i in Eq. (TT5T) . This strategy have been adopted in [11], where also all needed 
integrals have been carefully classified and computed. 



7 



Although rather transparent, the above derivation of Ri has two drawbacks. First of 
all, it requires the knowledge of the spurious terms [§. Secondly, it is not suitable when 
combining diagrams together because, when taking common denominators, additional 
terms containing q 2 appear in the numerator, that may give rise to new rational parts. 
The bookkeeping of such new structures is equivalent to the treatment of each diagram 
separately, jeopardizing the ability of the 4-dimensional OPP technique of dealing directly 
with the amplitude. For these reasons we present, in the next section, a different way of 
attacking this problem that does not relies on spurious terms and that also allows one 
to combine diagrams before fitting the coefficients d, c,b,a. This second method is better 
suited for a numerical implementation, and it has been already successfully implemented 
in a Fortran code [12]. 

4 The n-dimensional cuttings to compute Ri 

The Rational Terms Ri can be computed by looking at the implicit mass dependence 
(namely reconstructing powers of q 2 ) in the coefficients d, c, b of the one-loop functions, 
once q 2 is reintroduced through the mass shift 



This procedure is formally equivalent, in the generalized unitarity framework, to the appli- 
cations of n-dimensional cuts, and is obtained, in the OPP language, by simply performing 
the OPP expansion of Eq. (fTPl directly in terms of the n-dimensional denominators of 
Eq. f fTB . By doing so, all coefficients of the OPP expansion start depending on q 2 . The 
spurious terms keep being spurious, because they vanish due to Lorentz invariance (that is 
untouched when including powers of q 2 ), while the coefficients d, c, b generate the following 
extra integrals [10] 




(16) 



/ 



in 2 2 2 

~Y m i + m j 



3 



2" 



d n q 



DiDj 



+ 0(e) , 




(17) 



One can prove that 



b(ij;q 2 ) 
c(ijk; q 2 ) 



b(zj) + q 2 b^(zj), 
c(ijk) + <fc™(ijk) . 



(18) 



2 After multiplication with the Zj, they give non vanishing contributions. 
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Furthermore, by using Eq. (IT6j ). the first line of Eq. (ITDl becomes 

m—1 m—1 

v^ m \q,q 2 )= Yl [%i¥ 3 ;g 2 ) + 4¥^3;9 2 )] II A, (19) 

and the following expansion holds 

m 

pM ( g , f) = g (2j - 4) rf (2j '- 4) (?) , (20) 
i=2 

where the last coefficient is independent on q 

d^ 2m - 4 \q) = d (2m - 4 K (21) 

In practice, once the 4-dimensional coefficients have been determined, one simply redoes 
the fits for different values of g 2 , in order to determine b^ 2 \ij), 6 2 \ijk) and d <y2m ~ 4 \ Such 
three quantities are the coefficients of the three extra scalar integrals listed in Eq. (ITTj) . 
respectively, so that 

«0<H<*2 

-iE & (2) (Mi) (< + < - ^nl!) . (22) 

io<ii V / 

In Appendix [Al we prove Eqs. (IT8j)- (l2~Tj) , we single out the origin of d^ 2m ~^ as the coefficient 
of the last integral of Eq. ( ITTj) and we show how it can be also derived outside the OPP 
technique. Finally, yet another way of computing d^ 2m ~^ can be obtained by noticing 
that 

d C*n-*) _ lim PM M 2 ) (23 ) 

This last method is also implemented in the code of Ref. [12]. 

We close this section by stressing that the way of computing the coefficients appearing 
in Eq. ( 1221 is immaterial. Therefore the method to extract Ri described in this section, 
namely by looking at the mass dependence of the coefficients of the scalar loop functions, 
can be used independently on the OPP technique. In particular, one can derive all needed 
coefficients also with the help of analytical methods. 



5 Conclusions 

We have discussed and clarified the origin of the Rational Terms appearing in one-loop 
amplitudes, showing that they can be classified in two classes. The first class (R2) can be 
computed by defining tree-level like Feynman rules for the theory at hand. We precisely 
outlined the way to derive the needed extra Feynman rules, listing them explicitly in the 
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case of QED. We therefore consider the problem of computing R 2 completely solved. The 
second piece (Ri) can be calculated in two different ways. We presented a first technique 
that relies on the OPP method and a second, more general, computational strategy. Both 
methods have been successfully tested within the OPP method. The second one, however, 
is more suitable for a numerical implementation, and it has been used in the numerical 
code CutTools. 
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Appendices 

A The q 2 dependence of the OPP coefficients 

Our starting point is a rank r tensor m-point integrand defined as 

A m , r = ^ ' ' (m>3,r>l). (24) 

By expressing the integration momentum q in terms of the basis of the external vec- 
tors, with coefficients linearly depending on the propagator function appearing in the 
denominator [10], one ends up in an expression that it is identical to the OPP master 
equation plus terms containing higher-point scalar amplitudes A^.q k = 5, . . . ,m. Since 
all reductions formula used so far are simple polynomial in terms of q 2 , the same is true 
for all coefficients appearing in that expansion. This, in conjunction with the fact that 
the maximum rank allowed r — m, easily proves first of all Eq. ( 1T81 ). In the next step 
one reduces the k— point scalar terms k — 5, . . . ,m in terms of the 4— point ones at the 
integrand level. After that step the individual d and d coefficients become rational func- 
tions of q 2 . Nevertheless, since the £> (m) (g,g 2 ) of Eqs. (pi) and ((201) is nothing more than 
the combined numerator of all scalar terms with k = 4, . . . , m, and all coefficients are 
polynomials in q 2 , so is the V^ m \q,q 2 ) function. Finally it is straightforward to see that 
in the OPP expansion in terms of n— dimensional propagators, the q 2 — > oo behavior of 



10 



the individual d terms is g 4 , which means that the only integral involved is the last one 
of Eq. (HI). 

Another way to derive the same results is by using the reduction at the integrand level 
introduced in [13]. One can express A m . r as a linear combination, with tensor coefficients, 
of five classes of lower rank tensors: A m;r _i, q 2 A m . r _2, An ; r-2, An-ijr-i an d v4 m _i ;r _ 2 . To 
keep things as transparent as possible, we omit explicitly writing the coefficients and we 
denote such a linear combination using the following notation [f| 

A m ;r = {^m;! — \\q -A m ., — 2 |-*4m;j — 2 |-^m— l;i — 1 |-^m— l;r— 2 J • (25) 

Analogously, it was proven that 

An;l = {An;o\An-l;o} (m > 4) , 

A 4;1 = {A 4 .o|A3.oM(g)A 4; o} , (26) 
where d(q) is defined such as 

'd n qd(q)A 4 . = 0. (27) 



By using Eqs. (1251) and (1261) it is easy to constructively prove Eqs. (12011 and (T21T) . We 
explicitly give the derivation for the case m = 6 and r = 6. By iteratively applying 
Eqs. ( l25l ) and ( l26l) , one ends up with 

4i;6 = E {^;0|?%;0|<?%-;0} + {<f A;C)} + {%)^4;0 1 %)<f ^4;0 } 
.7=4,5,6 

+ 0(^3;r 3 ), (28) 

where 0(A 3 - r3 ) means that we are neglecting contributions with 3 or less denominators. 
By power counting, only the term g 4 A 4;0 contributes to R\. Notice also that its coefficient 
(that we call z 4 ) is independent on q. Now we can take a common denominator in Eq. (T28l) 
by multiplying and dividing 5 and 4-point structures by the relevant missing n-dimensional 
propagators. In particular, for example, by calling Di and Dj the 2 denominators that do 
not appear in y4 4;0 

z A q 4 A 4 . = z^DiDjAe-v = z A q\q 2 + A) (f + Dj) A 6;0 . (29) 

The numerator of the resulting expression is polynomial in q 2 and it is nothing but the 
function T>( m \q, q 2 ) of Eqs. (fT9l and (l20l . with m = 6. Furthermore d^ = z 4 , independent 
on q. The general case can be derived along the same lines. 

Eq. ( 1291) also clarifies why d^ 2m ~^ is the coefficient of the 4-point like last integrals 
of Eq. (1T7T) . In fact, m — 4 among the m original n-dimensional denominators always 



3 Notice that each of the five terms of Eq. I|25p may actually represent an entire class of contributions 
with different combinations of denominators. For example A m _i ;r _i stands for all m rank r — 1 tensor 
integrands that can be obtained by omitting 1 among the original m possible denominators. 
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completely factorize in front of d^ 2m ~^ . Notice also that the origin of the coefficient 
^(2m-4) - g un jq Ue iy coming from the mass dependence of the coefficients of the 4-point 
scalar functions after tensor reduction, but before expressing m-point scalar functions with 
m> 4 in terms of boxes. The reason why we do not reduce Eq. ( l28l) to structures with 4- 
denominators is that this would bring a q 2 dependence in the denominator, when passing 
from 5 to 4 denominators. In our notation [14] 



^5;0 — \ At;0 

I q 2 + Ci 



-d(q)A,. , (30) 



q 2 + Ci 

with Ci constants. It is therefore much better to take, instead, common denominators. 
Finally, analogous techniques can be used to prove Eq. (fT8l . 
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